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ABSTRACT 

We  present  three-dimensional  analysis  of  two-dimensional  guided  resonances  in  photonic 
crystal  slab  structures.  This  analysis  leads  to  a new  understanding  of  the  complex  spectral 
properties  of  such  systems.  Specifically,  we  calculate  the  dispersion  diagrams,  the  modal 
patterns,  and  transmission  and  reflection  spectra  of  these  resonances.  From  these  calculations,  a 
key  observation  emerges  involving  the  presence  of  two  temporal  pathways  for  transmission  and 
reflection  processes.  Using  this  insight,  we  introduce  a general  physical  model  that  explains  the 
essential  features  of  complex  spectral  properties.  Finally,  we  show  that  the  quality  factors  of 
these  resonances  are  strongly  influenced  by  the  symmetry  of  the  modes,  and  the  strength  of  the 
index  modulation. 

INTRODUCTION 

Photonic  crystal  slabs  are  a particularly  important  class  of  photonic  crystal  structures.  A 
photonic  crystal  slab  consists  of  a two-dimensionally  periodic  index  contrast  introduced  into  a 
high-index  guiding  layer  (inset  in  Figure  la).  These  structures  support  in-plane  guided  modes 
that  are  completely  confined  by  the  slab  without  any  coupling  to  external  radiations.  These 
guided  modes  allow  the  control  of  light  within  the  layer  at  the  wavelength  scale.  Therefore,  the 
slab  structure  may  provide  the  basic  substrate  for  large-scale  on-chip  integration  of  photonic 
components  and  circuits.  [1-8] 
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Figure  1.  The  band  structure  for  (a)  even  and  (b)  odd  modes  in  a photonic  crystal  slab.  The 
structure  of  the  slab  is  shown  in  the  inset  of  (a),  and  consists  of  a square  lattice  of  air  holes  with  a 
radius  of  0.2a  introduced  into  a high-index  dielectric  slab  with  a dielectric  constant  of  12  and  a 
thickness  of  0.5a.  Even  and  odd  modes  are  defined  with  respect  to  the  mirror  parallel  to  the  slab. 
The  gray  regions  are  the  continuum  of  radiation  modes.  Solid  lines  outside  the  gray  region  are 
guided  modes.  Solid  lines  within  the  gray  region  are  guided  resonances. 
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In  addition  to  in-planc  wave  guiding,  photonic  crystal  slabs  can  also  interact  with  external 
radiations  in  complex  and  interesting  ways.  Of  particular  importance  here  is  the  presence  of 
guided  resonances  in  the  structures.  [9-15],  Similar  to  the  guided  mode,  a guided  resonance  also 
has  its  electromagnetic  power  strongly  confined  within  the  slab.  Unlike  the  guided  mode, 
however,  the  resonance  can  couple  to  external  radiations.  Therefore,  guided  resonances  can 
provide  an  efficient  way  to  channel  light  from  within  the  slab  to  the  external  environment.  This 
property  has  been  exploited  in  the  designs  of  novel  photonic-crystal  based  light-emitting  diodes 
[1 1,  16],  lasers  [17,  18],  and  directional  output  couplers  [ 19J.  In  addition,  the  guided  resonances 
can  significantly  affect  the  transmission  and  reflection  of  externally  incident  light,  resulting  in 
complex  resonant  line  shapes  that  are  useful  in  filter  applications.  [9][20] 

The  purpose  of  this  paper  is  to  present  a novel  analysis  of  guided  resonances  in  photonic 
crystal  slabs.  Our  analysis  elucidates  a variety  of  complex  spectra  phenomena  associated  with 
these  resonances.  We  compute  the  dispersion  diagrams  and  the  eigen-field  distributions  of  these 
resonances  with  a plane-wave  band  structure  computation  method.  We  then  perform  finite- 
difference  time-domain  simulations  to  determine  the  transmission  and  reflection  spectra,  and  to 
visualize  in  real  time  the  interaction  between  the  resonances  and  the  incident  light.  Emerging 
from  these  simulations  is  a key  insight  that  involves  the  presence  of  two  temporal  pathways  in 
the  transmission  and  reflection  processes.  Based  upon  this  insight,  we  introduce  a general  and 
intuitive  theory,  which  uses  only  interference  and  energy  conservation  arguments,  to  explain  all 
the  complex  features  in  the  spectral  line  shapes.  Finally,  we  analyze  the  angular  and  the 
structural  dependence  of  the  guided  resonances,  and  we  show  the  wide  ranges  of  tunability  in 
quality  factors  for  these  resonances. 

BAND  STRUCTURE  OF  THE  GUIDED  RESONANCES 

Since  the  spectral  features  of  the  guided  resonances  in  a photonic  crystal  slab  will  turn 
out  to  depend  critically  on  their  modal  properties,  it  is  helpful  to  begin  our  discussion  with  a brief 
overview  of  the  band  structure  properties  of  these  resonances.  Throughout  this  paper,  our  model 
system  will  consist  of  a square  lattice  of  air  holes  introduced  into  a dielectric  slab  (inset  in  Figure 
1 ).  The  thickness  of  the  slab  is  0.5 a,  and  the  radius  of  the  holes  is  0.2 a,  where  a is  the  lattice 
constant.  The  dielectric  constant  of  slab  is  12,  which  roughly  corresponds  to  the  dielectric 
constant  of  Si  or  GaAs  at  optical  wavelengths.  For  such  a structure,  because  of  the  translational 
symmetries  within  the  plane  of  the  slabs,  the  physical  properties  of  the  slabs  can  be  described  by 
a band  diagram  that  relates  the  frequencies  of  all  the  three-dimensional  modes  to  the  in-plane 
wave  vectors  [ I ][2],  The  band  diagram  can  be  computed  by  a pre-conditioned  conjugate  gradient 
minimization  of  a Maxwell  operator  expanded  on  a plane  wave  basis  [21]. 

For  our  model  system  as  shown  in  the  inset  of  Figure  1 a,  the  band  diagram  for  the  even 
and  odd  modes  are  plotted  in  Figure  1 (a)  and  (b),  respectively.  Modes  below  the  light  line  are 
bona-fide  guided  modes  with  infinite  lifetime,  in  spite  of  the  large  index  contrast  introduced  by 
the  air  holes.  The  guided  modes  above  the  light  line,  on  the  other  hand,  cun  couple  to  radiation 
modes  and  possess  a finite  lifetime.  These  modes  therefore  become  guided  resonances.  They  are 
called  “guided”  since  they  are  closely  related  to  the  guided  mode  bands  in  a uniform  slab,  and 
should  therefore  retain  significant  portions  of  the  electromagnetic  power  within  the  dielectric 
slab. 

The  presence  of  the  air  holes  in  the  crystal  also  generates  a discrete  translational  and 
rotational  symmetry,  and  thereby  dictates  the  degeneracy  of  the  bands.  At  most  7-points,  (except 
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Figure  2.  The  frequencies  of  the  resonances  at  T as  a function  of  the  radius  of  the  holes  in 
the  slab.  The  slab  has  a dielectric  constant  of  12  and  a thickness  of  0.5a.  The  modes  are  four-fold 
degenerate  in  the  structure  without  holes.  For  structures  with  holes,  the  four-fold  degeneracy  is 
broken,  resulting  in  a pair  of  doubly  degenerate  states,  and  two  singly  degenerate  states. 


for  the  special  points  T,  X and  M),  the  bands  are  singly  degenerate.  At  the  T point,  the  point 
group  supports  a two-dimensional  irreducible  representation,  allowing  for  the  existence  of 
doubly  degenerate  states.  Therefore,  the  four-fold  degeneracy  at  the  T point  for  a uniform  slab 
splits  in  the  presence  of  the  air  holes,  as  clearly  seen  in  Figure  2,  where  we  plot  the  frequencies 
of  the  resonant  modes  at  F as  a function  of  the  radius  of  the  holes.  As  the  radius  of  the  holes 
increases,  the  modes  separate  into  a pair  of  doubly  degenerate  states  and  two  singly  degenerate 
states. 

For  the  crystal  structure  with  r = 0.20 a,  we  show  the  power  density  distribution  of  the 
first  resonant  band  at  T in  Figure  3.  The  mode  is  singly  degenerate  with  a frequency 
<y  = 0.35  • (2.71c /a).  Since  any  singly  degenerate  mode  should  belong  to  a one-dimensional 
irreducible  representation,  the  power  density  distribution  of  the  mode  should  possess  the  full 
symmetry  of  the  lattice.  This  can  be  seen  in  Figure  3(a),  which  shows  the  spatial  distribution  of 
the  power  density  on  a slice  parallel  to  the  slab.  Also,  the  resonant  nature  of  this  mode  is 
exhibited  in  Figure  3(b),  which  shows  that  the  power  density  is  strongly  confined  within  the  slab. 
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Figure  3.  Spatial  distribution  of  the  power  density  in  electric  fields  on  (a)  a horizontal  slice, 
and  (b)  a vertical  slice,  for  the  lowest-order  singly  degenerate  resonance  at  T.  The  lines  indicate 
the  position  of  the  interface  between  dielectric  and  air.  The  white  color  represents  low  intensity 
and  the  dark  color  represents  high  intensity,  as  indicated  by  the  color  bar  at  the  bottom  of  the 
figure. 
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The  band  structure  computations  thus  allow  us  to  examine  the  dispersion,  the  field 
distributions,  and  the  symmetry  properties  of  the  guided  resonances.  For  a complete 
understanding  of  these  resonances,  however,  we  must  also  study  their  lifetimes,  and  their 
interactions  with  external  radiations.  These  questions  will  be  addressed  in  the  next  session  by 
finite-difference  time-domain  simulations. 

TIME-DOMAIN  ANALYSIS  OF  THE  GUIDED  RESONANCES 

Computational  methods 

The  computational  domain  for  our  finite-difference  time-domain  study  [22]  includes  a 
single  unit  cell  of  the  crystal.  On  the  top  and  bottom  surfaces  of  the  computational  domain,  we 
impose  the  PML  absorbing  boundary  conditions  [23].  For  the  remaining  four  surfaces  that  arc 
perpendicular  to  the  slab,  we  impose  a Bloch  periodic  boundary  condition  on  the  electric  fields 

E: 


£(/+ a)  = /"£(/•)  (1) 

Here,  a is  a lattice  vector  of  the  square  lattice,  and  k is  a wavevector  that  is  parallel  to  the  slab. 
We  note  that  by  Bloch’s  theorem,  k is  a conserved  quantity  in  the  scattering  process. 

We  generate  an  incident  plane  wave  by  placing  a source  plane  consisting  of  oscillating 
dipoles  near  the  top  surface  of  the  computational  domain.  For  two  dipoles  in  the  plane  that  are 
separated  by  a distance  vector  r,  we  set  the  relative  phase  between  them  to  be  f . Therefore,  in 
combination  with  the  boundary  condition  as  specified  in  Eq.  ( 1 ),  the  source  plane  generates  an 
incident  plane  wave  with  a parallel  wavevector  component  k.  In  addition,  the  amplitudes  of  the 
dipole  moments  are  set  to  oscillate  at  a constant  frequency  with  a Gaussian  profile  to  create  a 
temporal  pulse.  This  computational  setup  thus  allows  us  to  calculate  the  response  functions  of 
the  structure  at  a given  k for  a wide  range  of  frequencies  in  a single  simulation  run.  (Notice  that 
this  is  not  a constant  incidence  angle  calculation.  At  a fixed  parallel  wavevector  k,  the  incidence 
angle  changes  with  frequency.) 

The  transmission  and  reflection  spectra  are  obtained  by  first  Fourier  transforming  the 
recorded  time  sequence  of  field  amplitudes  at  their  respective  monitor  points.  (We  note  that 
monitoring  the  field  amplitudes  only  at  the  two  monitor  points  is  valid  for  the  frequency  range 
1 0<2nc!a , where  no  diffraction  occurs.)  The  spectra  are  then  normalized  with  respect  to  the 
incident  pulse  that  is  calculated  in  an  identical  simulation  in  vacuum  without  the  slab  structure. 
For  reflection,  the  normalization  step  is  preceded  by  subtracting  the  incident  pulse. 

Transmission  and  reflection  spectra 


Using  the  computational  setup  as  described  in  Section  III.  1,  we  calculate  the  transmission 
and  reflection  coefficients  at  various  k-points  for  the  structure  as  shown  in  Figure  1 . In  the  case 
where  k = x ■ 0.2-  2/r la , the  calculated  spectra  for  the  .s-polarizcd  incident  wave  arc  shown  in 
Figure  4.  (An  .s-polarized  wave  has  its  electric  field  perpendicular  to  the  plane  of  incidence.  In 
this  case,  the  electric  field  is  polarized  along  the  y-direction).  The  spectra  consist  of  sharp 
resonant  features  superimposed  upon  a smoothly  varying  background. 
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Figure  4.  (a)  Transmission  and  (b)  reflection  spectra.  The  solid  lines  are  for  the  photonic 

crystal  structure  shown  in  the  inset  of  Figure  1(a).  The  dashed  lines  are  for  a uniform  dielectric 
slab  with  a frequency  dependent  dielectric  constant,  as  defined  in  Eq.  (2),  and  a thickness  of 
0.5a.  The  incident  wave  is  .v-polarized,  and  has  a parallel  wavevector  k = 0.2  ■ 2k  la  along  the  x- 
direction. 

The  background  in  Figure  4 resembles  the  Fabry-Perot  oscillations  when  light  interacts 
with  a uniform  dielectric  slab.  To  clearly  see  this,  we  fit  the  background  to  the  spectra  of  a 
uniform  slab,  which  are  shown  as  dashed  lines  in  Figure  4.  The  uniform  slab  has  the  same 
thickness  of  0.5a  as  the  crystal,  and  the  light  is  incident  with  the  same  polarization  at  the  same 
parallel  wavevector.  The  dielectric  constant  of  the  uniform  slab  £i,  as  obtained  by  the  fitting 
procedure,  represents  an  effective  dielectric  constant  for  the  photonic  crystal.  Due  to  the 
presence  of  the  holes,  such  £i  is  a slowly  varying  function  of  the  frequency.  At  low  frequencies, 
the  wavelength  of  incident  light  is  large,  and  £|  for  this  polarization  approaches  the  average 
dielectric  constant  £avsof  the  crystal.  At  higher  frequencies,  as  the  incident  wave  probes  more 
details  of  the  crystal  structure,  £i  starts  to  deviate  from  £avg.  Within  the  frequency  range  in  Figure 
4,  i.e.  between  0.25  • 2m:  la  and  0.60  ■ 2 kc  la,  we  have  found  that  a frequency-dependent 
dielectric  constant 

£;  (to)  =-14.16  - ®2  + 15.18  -01+7.18  (2) 

gives  a very  good  fit  of  the  background  (Figure  4).  The  fit  here  corresponds  to  varying  CD  from 
10.62  at  to= 0.25  • 2 kc  la  to  1 1.5  at  to =0.60-  Inc  la.  (As  a comparison,  the  average  dielectric 
constant  for  the  crystal  is  10.6.)  Therefore,  except  for  the  sharp  resonance  features,  the 
background  of  the  spectra  for  the  crystal  can  be  adequately  accounted  for,  using  the  model  of  a 
uniform  dielectric  slab  with  a frequency-dependent  dielectric  function. 

Line  shape  analysis 

Superimposed  upon  the  smooth  background  in  the  spectra  for  the  crystals  are  sharp 
resonant  features.  Such  features  come  from  the  guided  resonances  of  the  slab.  In  most  cases,  the 
line  shapes  for  these  resonances  are  asymmetric  and  rather  complicated.  Extensive  experimental 
and  theoretical  work  has  been  performed  for  guided  resonances  in  structures  with  one- 
dimensionally  periodic  index  variation.  [24-31]  For  structures  with  two-dimensional  periodicity, 
these  resonances  have  also  been  studied  numerically  using  the  Rigorous  Coupled  Wave  Analysis 
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(RCWA)  method  [20],  and  analytically  using  vector  coupled-mode  theory  [ 1 3] [ 1 5].  Here,  we 
would  like  to  present  a novel  analysis  from  a time-domain  perspective.  We  will  observe  the 
important  features  in  the  time-domain  signatures  of  the  resonances.  And,  based  upon  the 
observation,  we  will  introduce  a general  and  intuitive  model  to  account  for  the  underlying 
physics. 

The  transmission  and  reflection  spectra  are  related  to  the  time-varying  fields  by  a Fourier 
transformation.  It  is  therefore  informative  to  examine  the  time  dependency  of  the  fields.  As  an 
example,  we  show  in  Figure  5(a)  the  electric  field  amplitude  at  the  transmission  monitor  point  as 
a function  of  time  steps,  for  the  calculation  that  gives  the  spectra  shown  in  Figure  4.  The  time 
sequence  consists  of  two  distinct  stages:  an  initial  pulse,  and  a tail  of  long  decay. 

The  presence  of  these  two  stages  indicates  the  existence  of  two  pathways  in  the 
transmission  processes.  The  first  pathway  is  a direct  transmission  process,  where  a portion  of  the 
incident  energy  goes  straight  through  the  slab  and  generates  the  initial  pulse.  The  Fourier 
transformation  of  the  initial  pulse  should  account  for  the  background  in  the  transmission  spectra. 
The  second  pathway  is  an  indirect  transmission  process,  where  the  remaining  portion  of  the 
incident  energy  excites  the  guided  resonances.  The  power  in  the  resonances  then  decays  slowly 
out  of  the  structure  and  produces  the  long  decaying  tail.  By  Fourier  transforming  the  decaying 
tail,  we  obtain  the  typical  symmetric  Lorentzian  line  shapes,  as  shown  in  Figure  5(b).  The 
analysis  of  the  resonant  line  shape  thus  allows  us  to  determine  the  quality  factor  Q of  the 
resonance.  A few  examples  of  the  Q values  for  this  structure  are:  360  for  the  resonance  at 
ro=  0.3 1 ■ (2/rc/a),  and  2500  for  the  resonance  at  co  = 0.36  ■ (2 Jtc/a).  The  interference  between 
the  direct  and  the  indirect  pathways,  therefore,  determine  the  transmission  property.  The  same 
observation  can  be  made  for  the  reflected  amplitude  as  well.  In  solid  state  and  atomic  physics, 
similar  interference  phenomena  are  commonly  referred  to  as  the  Fano  resonances  [32],  Such 
temporal  interference  phenomenon  has  also  been  analyzed  previously  for  surface  plasmon  in 
metallic  thin  films  [33], 

Taking  into  consideration  the  interference  between  these  two  pathways,  we  can  construct 
a simple  and  intuitive  model  that  quantitatively  explains  the  line  shape.  We  express  the 
transmitted  amplitude  r,  and  the  reflected  amplitude  r,  as  follows: 


t = tj+f 


i((0-CO„)+f 


Here,  t,i  and  ra  are  the  direct  transmission  coefficients.  01,  and  y are  the  center  frequencies  and  the 
widths  of  the  Lorentzian  from  the  resonance,  and  the  factor/ is  the  complex  amplitude  of  the 
resonant  mode. 

The  plus/minus  sign  in  Eq.  (4)  corresponds  to  resonant  modes  that  are  even/odd  with 
respect  to  the  mirror  plane  parallel  to  the  slab.  We  note  that  the  Lorentzian  functions  in  Eqs.  (3) 
and  (4)  correspond  to  the  decaying  amplitudes  of  the  resonances  to  the  reflection  and 
transmission  sides  of  the  slab,  respectively.  For  an  even  mode,  the  decaying  amplitudes  to  the 
two  sides  of  the  slab  are  in  phase,  while  for  an  odd  mode  the  decaying  amplitudes  are  180-degree 
out  of  phase.  Thus,  the  signs  in  Eq.  (4)  are  different  for  modes  with  different  mirror-plane 
symmetry  properties. 
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Figure  5.  (a)  The  field  amplitude  at  the  monitor  point  as  a function  of  time  step,  for  the 

same  calculation  as  shown  in  Figure  4.  Notice  the  existence  of  two  separate  stages:  an  initial 
pulse,  and  a long  decaying  tail,  (b)  The  Fourier  transformation  of  the  amplitude  as  shown  in  (a) 
from  time  step  20,000  - i 00,000.  The  spectral  intensity  exhibits  Lorentzian  line  shapes. 

The  factor/ can  in  fact  be  determined  purely  by  energy  conservation  arguments.  We  note 

that: 


W2+M2=1-  (5) 

Moreover,  since  ra  and  td  are  the  transmission  and  reflection  coefficients  through  an  uniform  slab 
with  the  appropriate  effective  dielectric  constant,  we  should  have 


Constraints  (5)  and  (6)  together  uniquely  determine  the  factor/.  Plugging  Eqs  (3),  (4), 
and  (6)  into  Eq.  (5),  we  have,  for  any  co: 


-2t/f 


_2L 


(co-a)0y  + f 


= tfh  ±rd\ 


J(cO-CD0f  + r 


COS 


arg(/)-arg(trf±/,)-aco: 


■J(co-co0)2 + f 


, (7) 


which  can  only  be  satisfied  if: 

f = <t,,±r„)-  (8) 

It  is  interesting  to  note  here  that  the  factor  f is  independent  of  the  resonant  line  width  y. 

The  parameters  r(i  and  td  represent  the  background  of  the  spectra.  Therefore,  as  discussed 
earlier  in  Section  III.2,  such  parameters  can  be  determined  by  fitting  the  background  to  the 
response  spectra  of  a uniform  slab,  as: 
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(9) 


k -k 

i— -sin  (£_,/?) 

MA-  ; 

cos(A„  ft)  - / sin  (A  ) 


,.2  ,,.2  • 

4\-0*:l 

for  a plane  wave  with  parallel  wavevector  kx,  incident  from  vacuum  with  a dielectric  constant 
fj,  = 1 , through  a uniform  dielectric  slab  with  a thickness  h and  a dielectric  constant  e [34],  The 
parameters  k M and  k,\  in  Equations  (9)  and  (10)  represent  the  wavevector  components  along  the 
z-axis  in  the  uniform  slab,  and  are  defined  as: 


In  obtaining  Equations  (9)  and  (10),  we  assume  a positive  frequency  convention,  in  order  to  be 
consistent  with  the  Lorcntzian  functions  that  we  have  chosen  for  the  resonance  in  Equation  (3) 
and  (4). 

We  note,  in  particular,  when  r,  = 0,  and  /,  = 1 , from  Eqs.  (3),  (4),  and  (8).  the  reflection 
and  the  transmission  coefficients  become: 


ijeo-co,,) 
i((0- a>„)+  y 


and 


i{(o-(o„)+y 


(13) 


(14) 


The  line  shapes  thus  become  symmetric,  and  the  structure  behaves  as  a narrow-band  reflector 
with  a Lorentzian  reflectivity  line  shape.  Previously,  this  scenario  was  noted  by  Wang  and 
Magnusson  [26].  In  the  general  case  when  * 0,  on  the  other  hand,  the  line  shape  becomes 
asymmetric.  The  transmission  can  vary  from  0%  to  100%  within  a very  narrow  frequency  range. 
A small  shift  in  the  resonant  frequency  may  therefore  lead  to  a drastic  change  in  the  response 
function.  This  effect  may  be  exploited  in  the  design  of  optical  sensors  and  switches. 
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Figure  6.  Comparison  of  theory  and  simulations.  The  empty  circles  in  (a)  and  (b)  are 
numerical  results  taken  from  Figure  5(a),  which  corresponds  to  the  two  lowest-frequency 
resonances.  The  solid  lines  are  theoretical  predictions  from  Eqs.  (3),  (8),  (9)  and  (10).  The 
parameters  of  the  theory  for  the  two  resonances  are:  (a)  ft)0  =0.3076-  Inc/a, 
jz=  4. 1 91  • 1 O'4  ■ 27ir  /a ; (b)  O)0  =0.3601- Inc/a,  y=  7.2483 ■ 1 O'5  • 2m.-/a. 


We  compare  our  theoretical  predictions,  as  defined  by  Eqs.  (3),  (4),  (8),  (9)  and  (10),  to 
the  numerical  results  for  the  first  two  resonances  shown  in  Figure  5.  (Both  of  these  resonances 
are  even.)  The  frequency  to  and  the  width  y of  each  resonance  are  determined  from  the 
simulations.  The  only  fitting  parameter  here  is  the  effect  dielectric  constant  Ei,  which  we  take 
from  Equation  (2).  The  theoretical  results  thus  obtained  are  shown  as  solid  lines  in  Figure  6.  The 
theory  agrees  completely  with  the  numerical  simulations. 

Wavevector  dependency  of  the  resonances 

To  explore  the  wavevector  and  polarization  dependency  of  the  resonances,  we  performed 
calculations  at  different  values  of  £x,  for  an  incident  wave  that  is  either  s-  or p-  polarized.  (The  s- 
polarization  has  the  electric  field  perpendicular  to  the  plane  of  incidence,  while  the  p-polarization 
has  the  magnetic  field  perpendicular  to  the  plane  of  incidence).  We  determine  the  position  and 
the  width  of  the  resonances  by  Fourier  transforming  the  decaying  tail,  as  discussed  earlier  in 
section  III.3.  The  results  are  summarized  in  Figure  7(a),  where  we  show  the  frequencies  of  the 
resonances  as  a function  of  kx.  Incident  waves  with  different  polarizations  excite  different 
resonances,  since  the  two  polarizations  possess  different  symmetries  with  respect  to  yz-mirror 
plane. 

We  note  that,  in  Figure  7(a),  some  of  the  bands  do  not  continue  to  the  T point.  In  other 
words,  certain  resonances  at  T do  not  couple  to  either  polarization  of  the  incident  wave.  A closer 
examination  of  Figure  7(a)  reveals  that  all  these  uncoupled  resonances  are  singly  degenerate. 
Previously,  this  effect  was  observed  experimentally  by  Pacradoni  et  al  [14],  and  discussed 
theoretically  by  Paddon  and  Young  [13],  and  by  Ochiai  and  Sakoda  [25] 

To  further  explore  the  wavevector  dependency  of  the  resonance,  in  Figure  7(b)  we  plot 
the  quality  factors  of  the  resonances  as  a function  of  kx,  for  the  four  lowest  bands  in  Figure  7(a). 
For  bands  with  different  symmetry  properties,  the  behavior  of  the  quality  factors  is  very 
different.  The  Q factors  approach  a constant  as  k „ vanishes  for  the  modes  that  connect  to  the 
doubly  degenerate  states  at  T.  For  the  modes  that  connect  to  the  singly  degenerate  states,  on  the 
other  hand,  the  Q factors  of  the  modes  diverge.  The  calculation  clearly  demonstrates  that  the 
symmetry  of  the  modes  can  significantly  influence  the  photon  lifetime  of  the  resonances. 
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Figure  7.  (a)  The  frequencies  of  (he  resonances  as  a function  of  wavevcctor,  for  the 

structure  as  shown  in  the  inset  of  Figure  1 a,  as  determined  from  the  time-domain  simulations. 
The  solid  circles  correspond  to  the  resonances  that  arc  excited  by  the  p-polarizcd  incident  waves, 
the  empty  circles  correspond  to  the  resonances  that  are  excited  by  the  s-polarizcd  incident  waves. 
Notice  that  some  of  the  bands  do  not  continue  to  T,  indicates  the  existence  of  uncoupled  states  at 
T.  (b)  The  quality  factor  as  a function  of  wavevcctor,  for  the  four  lowest  bands  in  (a).  The  solid 
lines  correspond  to  modes  that  connect  to  the  doubly  degenerate  state  at  T point.  The  broken 
lines  correspond  to  modes  that  connect  to  the  singly  degenerate  stale  at  T point. 


Frequency  (2  nc/a)  Frequency  (2 rtc/a) 

Figure  8.  The  transmission  spectra  at  normal  incidence,  for  crystal  slab  structures  with  a 
radius  of  (a)  0.05a;  (b)  0. 10a;  (c)  0.15a,  and  (d)  0.20a.  All  the  structures  have  a thickness  of  0.5a, 
and  a dielectric  constant  of  12. 
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Radius  dependence  of  the  resonance. 


In  addition  to  symmetry  related  effects,  the  lifetime  of  the  resonances  is  also  strongly 
influenced  by  the  radius  of  the  holes.  At  the  limit  where  the  radius  of  the  holes  approaches  zero, 
the  Q factor  for  all  the  resonances  should  diverge,  since  the  resonances  asymptotically  become 
real  guided  modes.  To  demonstrate  this  effect,  we  plot  in  Figure  8 the  transmission  spectra  at 
normal  incidence  for  four  different  structures  with  the  radius  varying  from  0.05a  to  0.20a.  The 
spectral  feature  for  the  resonances  indeed  becomes  sharper  as  the  radius  becomes  smaller.  For 
the  lowest-order  resonances,  the  Q factor  varies  from  approximately  5000  at  r = 0.05a,  to  213  at 
r = 0.20a.  At  a larger  radius,  the  Q-factor  should  be  even  lower. 

The  tunability  of  the  quality  factor  with  respect  to  the  radius  of  the  holes  is  important  for 
LED  and  laser  applications.  For  photonic-crystal  resonant-cavity  LED  structures,  optimal 
efficiency  occurs  when  the  line  width  of  the  resonances  become  comparable  to  the  line  width  of 
the  emitter  [35].  On  the  other  hand,  for  a laser  structure,  a high-Q  resonance  is  typically  desirable 
for  threshold  reduction.  Therefore,  as  we  have  demonstrated  in  this  paper,  photonic  crystal  slab 
structures  are  very  versatile,  and  can  be  specifically  tailored  for  different  light  emitting 
applications. 

SUMMARY 

In  summary,  we  present  a three-dimensional  frequency  and  time-domain  analysis  of 
resonances  in  photonic  crystal  slab  structures.  These  resonances  are  strongly  confined  with  the 
dielectric  slab,  and  yet  at  the  same  time  are  coupled  to  radiation  modes.  For  external  light 
incident  upon  these  slabs,  the  transmission  and  reflection  spectra  are  strongly  modified  by  the 
presence  of  these  resonances.  The  line  shapes  exhibit  complex  asymmetric  characteristics.  We 
show  that  all  the  complexities  in  the  line  shapes  can  be  accounted  for  with  a simple  analytic 
model  describing  the  interference  between  direct  transmission  (or  reflection),  and  the  exponential 
decaying  amplitudes  of  the  resonances.  We  also  demonstrate  that  the  quality  factors  of  these 
resonances  are  strongly  influenced  by  the  symmetry  of  the  modes,  and  the  radius  of  the  holes. 
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